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We construct a systematic mean-field-improved coupling constant and quark loop expansion for corrections to 
the valence (quenched) approximation to vacuum expectation values in the lattice formulation of QCD. Terms in 
the expansion are evaluated by a combination of weak coupling perturbation theory and a Monte Carlo algorithm. 
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The valence (quenched) approximation to the 
infinite volume, continuum limit of lattice QCD 
gives values for hadron masses and for meson de- 
cay constants not far from experiment. Missing 
from these calculations, however, is an indepen- 
dent theoretical estimate of the error arising from 
the valence approximation. A possible method 
for determining the valence approximation's error 
is proposed in Ref. [|l|. In the present article, we 
describe an improved version of of the method of 
Ref. which we believe is likely to require less 
computer time. The expansion we describe can 
be adapted to any choice of quark action but will 
be given here only for Wilson quarks. 

In an exact treatment of QCD, virtual quark- 
antiquark pairs produced by a chromoelectric 
field reduce the field's intensity by a factor which 
depends both on the field's momentum and on 
its intensity. In the valence approximation this 
factor, analogous to a dielectric constant, is 
approximated by its zero-field-momentum zero- 
field-intensity limit 1^ . Our expression for the er- 
ror in the valence approximation to any vacuum 
expectation value consists of a sum indexed by 
a power of a mean-field-improved coupling con- 
stant 1^ and by a quark loop count. Each term 
in the expansion requires as input, in effect, the 
dielectric constant entering the valence approxi- 
mation. The dielectric constant, for convenience, 
we obtain analytically from mean-field-improved 
perturbation theory to second order in the cou- 
pling constant. A related calculation without 
mean- field improvement is described in Ref. [Q. 
The remaining work of evaluating each term in 



the error expansion is done by a Monte Carlo al- 
gorithm. 

In a simple test case, our error estimate re- 
quires significantly less computer time than a di- 
rect comparison of full QCD and the valence ap- 
proximation. Whether this gain holds also for 
more interesting cases we do not yet know. 

We consider Wilson's formulation of euclidean 
QCD on some finite lattice, with periodic bound- 
ary conditions for gauge links u(x,y), and an- 
tiperiodic boundary conditions for the coupling 
matrix AI for a single quark flavor. For nj, either 
even or odd, degenerate flavors of quarks and any 
function of the gauge fields G, the vacuum expec- 
tation value obtained after integrating out quark 
fields is 



<G> 



z- 



dvdetiMYf expi-P), (1) 
6 

where, as usual, P is Wilson's plaquette action, fi 
is 6/(7^ for bare gauge coupling constant g, and v 
is the product of one copy of SU (3) Haar measure 
for each link variable on the lattice. We ignore for 
simplicity vacuum expectation values of products 
of quark and antiquark fields. 

The valence approximation to < G > of Eq. (|^) 



< G 

Z.., 



Z-' J duGexpi^P), 



dv exp(yP), 



(2) 
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with valence approximation given by Q/ gv^ for 
valence approximation bare gauge coupling g^ . It 
is convenient to define A/3 as the shift (3y — (3. 
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We now derive a coupling constant and quark 
loop expansion for the difference between a full 
QCD vacuum expectation value < G > and its 
valence approximation < G >i,. For this purpose, 
we recast < G > of Eq. (|^) as 



constant 



< G > 



z- 



dv G exp{^P - 



Q) 



Z = J di^expi^P + Q), 
Q = UftrlogiM) - ^P, 



(3) 



then introduce a parameter A multiplying Q, ex- 
pand < G > in powers of A, and replace A by 1. 
We obtain 

< G > = < G >,„ + (4) 

< (G- < G >v)iQ- < Q >„) >^ + . . . 

The nth term in Eq. (^) is the n — 1 quark 
loop correction to the valence approximation. As 
an error estimate for the valence approximation, 
we will concentrate on the one-loop correction 

< (G- < G >y){Q- < Q >v) >v. 

For tr\og{M) in the one-loop term of Eq. ^, 
we construct a coupling constant expansion. To 
do this, transform each gauge configuration to 
the euclidean lattice version of Landau gauge by 
maximizing X]y ^^['"(^i 2/)] each lattice site x. 
One algorithm for finding this transformation is 
discussed in Ref. We obtain some random 
sample of Gribov copies the details of which play 
no explicit role in our expansion. For each fixed 
gauge configuration, let z be the average over all 
lattice links of tr[u{x,y)]/3. Let Mq be a free 
hopping matrix with hopping constant kq chosen 
to give a quark mass 1/(2ko) — 4 that agrees with 
a mean- field- improved value of the quark mass 
1/(2zk) — l/(2zKc) carried by M, where k and Kc 
are, respectively, the hopping constant of M and 
the valence approximation to the critical value of 
this hopping constant. 

We now express trlog{M) in the form 



tr\og{M) = 

trlog{zMo[l-Mo'{Mo 



-'M)]}, (5) 



and expand to obtain a sum indexed, in effect, by 
powers of a mean-ficld-improvcd quark coupling 



Mog(M) = Mog(zA/o)- 
Y^^tr{[M~\Mo^ 



(6) 



-iM)]"}. 



The trace required in the second term of Eq. (^) 
we obtain from an ensemble of complex-valued 
quark fields (j)^ , 1 < r < R with each compo- 
nent of each field an independent random com- 
plex number on the unit circle. We obtain 

Mog(M) = trlog(zA/o)- (7) 
^-^ n 

nr 

where < . . . | . . . > is the inner product on the 
space of complex-valued quark fields. The term 
trlog(zMo) in Eq. (0) can be computed by diag- 
onalizing Mq in momentum space, and multipli- 
cation by Mq^ can be done efhciently using fast 
Fourier transformations. 

Combining Eqs. (||) and (0) gives a coupling 
constant and quark-loop expansion for correc- 
tions to valence approximation vacuum expecta- 
tion values. Still to be determined is the shift 
A/3. Eqs. (|^) and (Q) are formally correct for any 
choice of A/3. The rate at which these expansions 
converge, however, will be affected by the choice 
of A/3. One possible way to choose A/3 is to re- 
quire the error or some approximation to the error 
in a particular physical quantity to be zero. For 
the present discussion, we choose A/3 by requiring 
that (analytic) mean-field-improved perturbation 
theory give zero valence approximation error, to 
first order in quark loops and second order in the 
coupling constant, for the Landau gauge gluon 
propagator at minimal nonzero momentum. We 
obtain 



A/3 = 



9n/ 



4sin^(7r/L) < trll >„ 
[n22(p)-n22(0)] (8) 

n^-b) = -^^ir[r^(q+p/2)5(q+p) X 

T,{q+p/2)S{q)], 

where L is the lattice period, p is a momentum 
vector with a single nonzero component pi of 
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27r/L, < trll >„ is the valence approximation 
plaquette expectation value, and each component 
ofq in the sum over q ranges from tt/L to 27r— tt/L 
in steps of 2tt/L. The propagator S{q) has inverse 
S{qy^ of ^/i2Ko)-iY.^,l^lSin{qf,)-Y.^cos{qfJ,), 
and the vertex Tfj^(q) is siniji^) — ^cos{c[^) . The 
limiting value of A/3 for large L without mean- 
field improvement has been derived in Ref. [Q. 

As a test of our method we compared valence 
approximation expectations < G Eq. (||), 
their one-loop errors < (G— < G >i,)(Q— < 
Q >v) >v, Eq. (^), and the corresponding full 
QCD expectations < G >, Eq. (|^), for a lat- 
tice 10^ with (3y of 5.679, k of 0.16 and of 
2. According to Eq. (||), A/3 is then 0.243 giving 
a full QCD f3 of 5.436. For this case A/3 found 
by the method of Ref. is 0.244(6). We used 
224 independent equilibrium gauge configuration 
in the valence approximation ensemble, generated 
by an over-relaxed pseudo heat bath algorithm, 
600 random fermion fields R to evaluate the trace 
in Eq. (|^) and 176 weakly correlated equilibrium 
gauge configurations for the full QCD ensemble, 
generated by a red-black preconditioned hybrid 
Monte Carlo algorithm. The expansion in Eq. (|^) 
was carried to order n of 10. The calculation of 
< (G— < G >v){Q— < Q >v) >v was not turned 
carefully. In particular R of 600 in Eq. is much 
larger than its optimal vaule. The time required 
for the valence approximation and error calcula- 
tion was still less than 5% of the time required by 
the full QCD calculations. 

For G we used Wilson loops Wq , . . . Wio 
consisting, respectively, of paths 1x1, 2 x 
1, all rotations of steps in the directions 
1, 2, 3, —1, —2, —3, all rotations of steps in the di- 
rections i, 2, 3, -2, -i, -3, 3 X 1, 2 X 2, 4 X 1, 5 X 1, 
3 X 2, 4 X 2, and 3x3. For all loops, we found 
that the predicted error < (W— < Wi >v){Q— < 
Q >v) >v converges as a function of the highest 
power n of coupling strength used in Eq. (^) by 
n of 7 or smaller. For n of 7, Fig. |l] shows the 
relative shift of the valence approximation from 
full QCD {< > - < W^ >„)/ < W^ > and 
the predicted value < {Wi— < Wi >v){Q~ < 
Q >v) >v / < Wi >. To within statistical uncer- 
tainties, the predicted errors agree with the true 
errors. 
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Figure 1 . The predicted relative shift of 1 1 Wilson 
loops from their valence approximation values in 
comparison to the true shift of full QCD. 

The true errors in Fig. ^ were found from the 
shortest full QCD run sufficient to confirm equi- 
hbration of < Wq >,... < Wio >■ Nonetheless 
the statistical uncertainties in the predicted er- 
rors are much larger than those in the true er- 
rors. If we were to run the error prediction al- 
gorithm long enough to obtain statistical uncer- 
tainties comparable to the uncertainties found by 
a direct comparison of full QCD and the valence 
approximation, it is possible that the computer 
time required by the error algorithm would be- 
come comparable to that for full QCD. To find 
the uncertainty arising from use of the valence ap- 
proximation, however, the statistical uncertainty 
in the error estimate does not need to be too much 
smaller than the error estimate's central value. 
Used in this way, for the set of parameters of the 
test, our algorithm takes significantly less time 
than the shortest possible direct comparison of 
the valence approximation and full QCD. 
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